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Abstract 

We construct an Imbimbo-Mukhi type matrix model, which reproduces exactly 
the partition function of CP 1 topological strings in the small phase space, Nekrasov's 
instanton counting in N = 2 gauge theory and the large N limit of the partition 
function in 2-dimensional Yang-Mills theory on a sphere. In addition, we propose a 
dual Stieltjes-Wigert type matrix model, which emerges when all-genus topological 
string amplitudes on certain simple toric Calabi-Yau manifolds are compared with 
the Imbimbo-Mukhi type model. 
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1 Introduction 



Recently, classical geometries prove to be very powerful in determining strong coupling 
gauge theory dynamics. A familiar story is the AdS/CFT correspondence in M = 4 
SYM[T|. Similarly, computing the confining phase superpotential in M = 1 gauge theory 
has been facilitated by a large N duality for conifolds, called the geometric transition.]. 
On the dual deformed conifold, period integrals over holomorphic cycles are responsible 
for the effective superpotential of gluino condensates [3], HI E]. In this N = 1 context, 
Dijkgraaf and Vafa conjectured that the all-genus free energy of certain target CY's, 
locally diffeomorphic to T*S 3 , reduces to a hermitian (DV) matrix model[6j [7J [8]. As a 
further development of this celebrated matrix/geometry correspondence, for more general 
non-compact CY's, Aganagic et al. provided an unifying way in [9] , where the W-algebra 
symmetry (holomorphic diffeomorphism of CY's) in disguise is fully exploited to solve the 
all-genus free energy. In terms of the matrix model, these are known as W-constraints 
on the r-function of integrable hierarchies. Adopting the viewpoint that the partition 
function on CY 3-folds can be regarded as a multi-point function of non-compact B-branes, 
we can find its connection to the matrix model r-function via the Frobenius-Kontsevich- 
Miwa transformation. 



'unction of Af = 2 
type matrix model 



Our situation here is that the full non-perturbative partition 
Seiberg-Witten theory§| is found to be captured by an Imbimbo-Mukh: 
in the 't Hooft limit. Alternatively, because this instanton partition function |1 11 [T2l [T3] 
is precisely that of CP 1 topological strings in the small phase space (see, e.g. [HI 115]). 
our proposal may serve as another sort of matrix/geometry correspondence, though now 
the integrable structure is quite simple. 

The outline of this article is as follows. In the next section, we review the derivation 
of the Gromov-Witten prepotential of CP 1 topological strings in the small phase space. 
We focus on its free fermion representation. Though the integrable structure is simple, it 
is still fruitful enough to exhibit the underlying Toda hierarchy. In sec. 3, we construct an 
Imbimbo-Mukhi type matrix model to compute the CP 1 partition function and sketch a 
proof using essentially properties of Schur polynomials. Meanwhile, our result gives spon- 
taneously Nekrasov's instanton counting in M = 2 gauge theory. Comments on relations 
to 2-dimensional Yang-Mills theory are summarized in sec. 4. Based on observations of 
all-genus amplitudes of topological strings on some simple toric CY's, we also propose a 
Stieltjes-Wigert type matrix model dual to the previous Imbimbo-Mukhi type one in the 
't Hooft limit. 



1 See [2] for a recent review. 

2 It can be either abelian or non-abelian. The former case is realized by wrapping a D5-brane on 
a CP 1 C K3, where the non-perturbative effect arises from worldsheet instantons. By rearranging the 
abelian partition function, one as well obtains a non-abelian version, see Appendix B. 

3 The Imbimbo-Mukhi matrix model[10 is a r-function satisfying the 2-time Toda lattice hierarchy. 



1 



2 r-function of CP 1 model 

The complex projective space CP 1 is the simplest target manifold with a positive first 
Chern class. Since CP 1 has only two de Rham classes, i.e. the identity and the Kahler 
class, the topological A-model on it will have two corresponding physical observables, 
i.e. P and Q, respectively. Note that t Hi p and t n> Q (n > 0) indicating couplings to 
gravitational descendants of P and Q are identified with times of the Toda hierarchy in 
[T5| [T6] . When all couplings are turned off (in the small phase space) , the (genus zero) 
free energy is 

^o,p*o,Q + e*°^ (2.1) 

where £o,Q := t denotes the Kahler parameter of CP 1 and £q,p will pl a Y the role of the 
complex scalar vev in the J\f—2 U(l) vector multiplet. Besides, the exponential term 
comes from the degree-1 instanton. 

Following O [T7] , the partition function (r-function) can be expressed as a sum over 
Young tableaux A = (ki, ■ ■ • ,kj,---) (ki : i-th row length, |A| : number of boxes): 



z cp- _ V m|_ / ch 2 (a, A) 



(2.2) 



mA = n htaj = p = n h )H l > m ■■ ^ length. 

Here, h denotes the genus expansion parameter. Note that m A and d\, respectively, are the 
Plancherel measure and the dimension of A as a representation of S\ (symmetric group). 
In addition, the Chern polynomial Ch n is defined as a Fourier transform of the second 
derivative of f\(x, a) , which is the profile function obtained by rotating counterclockwise 
the Young tableau by |7r: 

1 r 00 u n 

-jdx f>(x,a) e u * = £ — ch„(a, A), (2.3) 



n=0 



where (box size: h x h) 

fx(x) = \x-a\+ p(x) 



p(x) = ^2 (\x — a — hki + h(i — 1)| — \x — a — hki + Hi\ m 4) 



+ \x — a + hi\ — \x — a + H(i — V 

Alternatively, Z cpl can be written into a more compact matrix element form: 

Z cpl = (0|e£ e tLo e^lO). (2.5) 
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To explain the notation, let us make use of fermionic operators 6, c satisfying the usual 
properties: 



{b r ,c s } = 5 r+s , , r,seZ + l/2; 6 r |0) = c s \0) = 0, Vr, s > (2.6) 

to define an orthogonal basis {|A)} w.r.t. A: 

oo 

|A) = IIW|I '). (A|«> = ^; c r |0') = 0, Vr, (2.7) 

i=i 

which satisfies 

e ^ 0) ^mitm w (2 - 8) 

and 

L |A) = ^ch 2 (a,A)|A). (2.9) 

Note that the Virasoro generator L = | J| + J-\J\ consists of the £7(1) current defined 
using the bosonization: 

J(z) =: b(z)c(z) := ^2 J n Z 1 , J n = '■ b s C n -s [Jmi J-n) = n ^m+n,0- 

n&L S&Z+1/2 

(2.10) 



For J = |, the matrix element reads 

= PYn , . 

h 2 \ 2 

It can be verified that (h = 1) 



Z cpl = exp±-(e t + -ta 2 ). (2.11) 



dt (Z CP (a)) 

which is none other than the first equation of the Toda lattice hierarchy. 



3 Imbimbo-Mukhi type matrix model 

Before writing down our matrix model, let us first study an illuminating generalization, 
namely, the Imbimbo-Mukhi matrix model which describes the deformed conifold T*S 3 
with B-brane insertions. 

According to [9], on the B-model side, to vary the complex structure of CY's is done 
by wrapping non-compact B-branes on degenerating fibers. The problem hence boils 
down to dealing with the moduli space of a Riemann surface. Moreover, with the special 
geometry relation, quantizing the complex moduli is realized by introducing quantum 
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Kodair a- Spencer fields. The vacuum expectation value of the KS field in turn defines a 
quantum Riemann surface which determines the CY with varied complex moduli. For 
full details on this subject, readers are recommended to consult [9]. 
We quickly review the case of T*S 3 , whose defining equation is 

uv — H(x,y) = 0, H(x,y) — xy — (3.1) 

Two KS fields (chiral bosons) 

x = d y (p(y), y = -d x (f){x) (3.2) 

are defined on asymptotic regions x, y — > oo of the degenerating locus H (x, y) = which 
has a topology of the sphere. If there are initially N non-compact B-branes at rc-patch 

(N\1>( Xl ).--i/>(x N )\W), ^ = eTs^\ (3.3) 

where \W) G 7i (TC: Hilbert space of chiral bosons 0) and (N\ denotes an iV-fermion 
state, the complex structure of T*S 3 is deformed to be 

a N 

xy = a _^ x{y + Y. nt nX n ) =», *n = --X>f ft . (3-4) 

n>0 71 i=l 

Notice that the canonical transformation relating x-patch and |/-patch gives rise to the 
iS-matrix, which acts on the B-brane operator as a Fourier transform, i.e. 

~ I f x y 

ip{y) = Sip(y) = / dx e~~s^ip{x). (3.5) 



In the deformed geometry (13.41) . moving B-branes from x-patch to |/-patch Y = diag(yi, • • • , y^) 
leads to 

f Y[dxi e-Ef^^M^^)...^^)). (3.6) 

^ i=l 

Further, via the Harish-Chandra-Itzykson-Zuber formula 

/ DU e iTv(BX) oc I , A . x , A : VanderMonde, (3.7) 

J A(x)A(p) 

where X = diag(xi, • • ■ ,xn) and B = U'PU with P = diag(pi, • ■ ■ ,Pn), we are able to 
obtain a Kontsevich-like matrix model for (13 .61) . i.e. 

Z(t,t) = — ^ / \\dxi A(x)(detXy^e^ ( ^ tnTrXn+ ^ x ^ ] 
A[y) J xx 

(3.8) 



JNxN 



where t n = ^-TrF n . For large enough N, we choose \i = g s N. (13.81) is the so-called 
Imbimbo-Mukhi matrix model first derived in 110!. 
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3.1 Proof 

We are now in a position to propose our model. Under the 't Hooft limit: 

fj, = Ng s : large but fixed, iV — > oo, g s — ► 0, (3.9) 

the following hermitian matrix integral reproduces Nekrasov's instanton counting formula 
as well as the previous Z cpl , i.e. 

Z(U) = / dX(detX)-^expTrJ V K n A n + XIT 1 Y (3.10) 



n>0 



where 



/iK n = --TrU n , U = diag(gV-. g = e _9s . (3.11) 

n 

Let us now sketch a proof of our proposal. A fact of particular use is that the r- 
function (I3.8P admits a S-matrix formulation [18], where the S-matrix has been explicitly 
diagonalized thanks to tachyon reflection amplitudes in the bosonic c = 1 string. Based 
on those, we rewrite Z(U) into (S: S'-matrix) 

00 TrU n 

Z(U) = (W>, \£) = expY. «-n|0), £ n = TrV n . (3.12) 

Given all mathematical details in Appendix A, we insert the identity operator 1 = 
Er\R){R\ into fl3TT2|) to get 

Z(U) =J2Mrs r (£)s r (£), Mr = Y[(ji-i+j), (3.13) 

R U(i,j) 

where Sr(£) is the Schur polynomial and the diagonal element Mr can be found in |19^l . 
Upon imposing the 't Hooft limit, Z(U) reduces to a MacMahon function|26[ l2?pl : 

N=oo i 

z(u) = = n 7j — ^- ( 3 - 14 ) 

n=l (1-W n )" 

We have taken advantage of the Cauchy's formula 

E?A TrB = n r J -V' ( 3 ' 15 ) 

where a, (6j) denotes the eigenvalue of A (B). Another way defining the Schur polynomial 
in terms of q is 

8 R (£) = TrV = q-& M »*M JJ /2 \ 2 , ^,j): hook length. 

(*,i)eA y y 

(3.16) 



4 Note that the diagonal element Mr also appears in the context of 1/2 BPS correlators in TV = 4 
SYM, see H3IH]. 

5 The subleading term ^^■/il- R l~ 1 in Mr, where 2 J^p^ ^ {—i + j) = k(R), is omitted when compared 
to the leading one. 
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By expanding q in (I3.16P around small g s , together with (12.21) . it is seen that 

£^**0 a - E(4f m 2 ,. (3.17) 

R R, "s 

Amazingly, the partition function of CP 1 in the small phase space is reproduced upon 

g s -> h, fi^e l + -aH. (3.18) 

It is fine that a can be turned off to zero because it stands for the complex scalar vev in 
the J\f=2 U(l) vector multiplet below. 

3.2 Nekrasov's formula 

As is shown, our Imbimbo-Mukhi type matrix model proves to give 

E(£f»5 (3-19) 



under the 't Hooft limit. As a matter of 
function in Af = 2 U(l) gauge theory [T%J® 



act, (I3.19P is nothing but Nekrasov's partition 
Z mst (a,A,h) = Y,exp(-S[f,X]), 



A 



= - [ [ dx l dx 2 fx{xi)fx{x 2 )lh{x l -x 2 ,A), 

/ J Jxi>X2 



(3.20) 



where 



7b(x,A) : = ^~ 
ds 



A s r°° dt e~ tx 



s=0i 



F(s) Jo t 1 - 3 (e ht - l)(e~ nt - 1) 



(3.21) 



r(s) = / dt e~* 1f~ x . 
Jo 

In addition, jn( x , A) has a genus expansion w.r.t. h: 

oo 

7 ,(x,A) = 5:^- 2 79 (x,A) 

9=0 



(3.22) 



where an identity 



1 d 1 1 ^ B 2g 



(e*-l)(e-*-l) tfte*-l t 2 ^2g(2g-2 



+ E 5 -7^' 2 '- 2 (3.23) 



3 How to convert (|3.19[) into a non-abelian instanton partition function is summarized in Appendix B. 
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and the definition of the Bernoulli number 

r 00 B 

\ - J->n n 



— =5^" (3 - 24) 



e x n=0 



have been used. Through the substitution 

fj, = L 2 logA + A, (3.25) 

we yield [14J 

^(a,A^) = E(|)' A| mi (3.26) 
Notice again that choosing a = does not lose any generality. 



4 Comments 

So far, we have constructed an Imbimbo-Mukhi type matrix model. In the 't Hooft limit, 
it reproduces the partition function of CP 1 topological strings in the small phase space. 
Besides, it is none other than Nekrasov's instanton counting formula in this limit. 

Let us include two more applications. The first is that our model can recover the 
partition function of 2-dimensional Yang-Mills theory on a sphere. Another aspect is that 
the all-genus free energy of special toric Calabi-Yaus can be derived from our model. To 
this end, we further propose a dual Stieltjes-Wigert type matrix model. 

4.1 2-dimensional Yang-Mills 

The partition function of 2d U(N) YM on a genus g and area A orientable manifold 

isplES] 

Z 2d = Y^M^exp ( - ^C a (J2)), (4.1) 

where A is the 't Hooft coupling of 2d YM, C^i?) is the quadratic Casimir of the repre- 
sentation R (conventions inherited from (12. 2ft and [38] : the box size is H x h) 

n r N 1 

C 2 (R) = (Nkt + h{h - 2i + 1)) = J dxf^x) (—x 2 + —x 3 ) (4.2) 

and 

dimR = rj h ~ k J±A ~ 1 

l<i<j<N 3 1 

= exp ( - 1 j J dxdyf£(x)f£(y) lh (x -v) + \j dxf'^x) ( lh (x + KN) - lh {KN) 

(4.3) 
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Here, we have set a = and A = 1 for Jr(x, a) and 7/1(2, A). Note that the second term 
in the above second line serves as a cut-off which will be dropped out when N is taken to 
infinity. 

At the first sight, when N goes to infinity, dimR — > itir and C^i?) ~ N\R\, so the 
genus zero Z 2d is exactly of the form (13. f 91) . Let us check this point precisely. Combining 
everything, we yield 

Z 2d = exp ( - 29 f f dxdyf R {x)f R {y) 7s(x - y) 



r v ^ J Jx> y 

+ J dxfUx) (ln(x + UN) - lh (hN)) - M J dxf'^x) { £- x * + _L X 

(4.4) 

For large N, we turn out to get 



Z 2d -> 2 exp ( - / jf dxdyf£(x)f£(y) lh {x - y) - ^ / dx/£(s) x 2 ) . 



(4.5) 



Utilizing another equality for the profile function /p(x)[l2]: 



4ft 2 

in the genus zero case, we are left with 



R \ = TrJ dx ?k{x) x 2 , (4.6) 



ii 



In other words, we may have recovered the 2d YM partition function in large A" limit 
using our matrix model upon identifying^ jj = e~ XA l 2 . 

4.2 Toric Calabi-Yau's and unknot Wilson loops 

It is well-known that U(N) level t Chern-Simons gauge theory on S 3 is large A" dual to the 
topological A- model on a resolved conifold O(-l) + 0(— 1) — > P 1 . That is, the partition 
function on the resolved conifold (Kahler parameter t = g s N, q = e~ 9s ) 

. , e ~ nt 2m 
Z P i = exp(-J^), T = — ^, g, 

n>0 n[2sia(ng s /2) 



£ + N (4.8) 

is that of Chern-Simons under the 't Hooft expansion |24[ [25] . This is an open/closed 
duality realized via the geometric transition of conifolds. 



7 Note that both /1 and e^ XA l 2 can be complex. For /x, this is true in the context of the CY manifold. 
In the case of — XA/2, it plays the role of the gauge coupling of a 4-dimensional theory (e.g. D5-branes 
compactified on the S 2 of a resolved conifold) and can be complexified by turning on the theta angle. 
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As another example |26j. let us consider two U(N) Chern-Simons on two S 3 's, separated 
from each other by a complexified Kahler parameter t'. The Ooguri-Vafa operator [28] 
associated with two holonomies in this setup is 

Z(U, V, t') = exp ( ~e nt ' Tr U n Tr V n ) = £ e' 1 '^ Tr U Tr V (4.9) 

n>0 71 R R R 

with the normalized vev 

(Z(U,V,t')) = Y J e~ t '^W R {lC)W R {lQ ^Ml fl l Sfl ( Si = q^)\ 



R R, 
W R {K) = 8 R (xi = g^" 2 ^ 1 )/ 2 ), M = e~U N - 



(4.10) 



Note that W R ()C) stands for an unknot Wilson loop in S 3 . After the geometric transition, 
a smooth toric CY emerges which is characterized by three Kahler parameters, namely, 
t 1= t 2 = g s N and t = t' - (*^). The free energy - \og(Z{U, V, t')) is 



rt>0 7i ( 2sm(ng s /2) 

It is found that (14. 10p readily agrees with our result in (13. 17j) if q — > q^ 1 and M — > fi are 
carried out in (14.101) . 

4.3 Stieltjes-Wigert type matrix model 

Inspired by ( 14.111) . we are capable of proposing a Stieltjes-Wigert type matrix model which 
is dual to ( 13.101) . That is, 



-nti 



+ e~ nt2 + e~ nt (l - e~ ntl )(l - e~ n * 2 ) 



[ dFexpTr ( - -L(logF) 2 + l g(Y ® l mxm - l NxN ® V) 
JnxN \ Lq* 



vo\U(N) Jnxn V 2g. 

V = diag(ui 



(4.12) 



where Y is hermitian and the measure gives the usual VanderMonde. Without the second 
logarithm term inside the potential, (14.121) is orig inally designated [2H1 EDI EB E^fl to 
rewrite the partition function of U(N) Chern-Simons theory on S* 3 . 

To prove the equivalence to (I3.10p . we note that (14. 12[) dictates a B-model calculation 
of the non-compact B-brane correlator [35l [36] with their moduli V. Put it differently, 
(14.121) is a mirror description of the A-model resolved conifold, on whose external leg of 
the toric diagram a stack of probe Lagrangian branes is inserted. Besides, each component 



See [33l|34] for more details about Chern-Simons theory and Stieltjes-Wigert polynomials. 
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of V is chosen to be Vi = exp —g s {t — i + m + which immediately implies that these 
Lagrangian branes are all located at the same place and, therefore, trigger a geometric 
transition such that a new toric CY appears, referred to as a bubbling CY in |37j . 

Translating everything into the A-model language, one can utilize the melting crys- 
tal method to write down the Lagrangian brane scattering amplitude via the quantum 
dilogarithm L(v, q) as (Kahler parameter of P 1 : t = (N + m)g s ) 

Z crystal = Z pl (t)M(q) n(l - ft T^Ty 

i<j V 3 i=l H V i e > Q) 

(4.13) 

L(v, q) = II (1 - n k ) = exp ( £ ^- j , [k] = q k ' 2 - q~ k '\ 

which is essentially the same as f)4.12p up to an irrelevant overall factor and a framing 
choice|36j. Moreover, by doing so and inserting the moduli variable V, it is seen that 
— log ( Zcrystai/ M(q) 2 J coincides with (14. lip , except now [37] 

ti = g s m, t 2 = g s (N - m), t = gj. (4.14) 

In other words, given the result discussed in sec. 4. 2, the Stieltjes-Wigert type matrix 
model in (14.121) can be regarded to be equivalent to (I3.10P (up to MacMahon's) in the 't 
Hooft limit through m = y = N' — * 00, q — > q^ 1 and e* — > fj,. 

Acknowledgements 

The author would like to thank Feng-Li Lin, Yutaka Matsuo, Kazuhiro Sakai, Yuji Sug- 
awara, Ryo Suzuki, Masato Taki, Satoshi Yamaguchi for helpful discussions and, espe- 
cially, Tohru Eguchi for reading carefully the manuscript. Also, he gratefully acknowledges 
the hospitality of the string group of National Taiwan University, where the final stage of 
this work was done. This work is partly supported by the 21COE program of University 
of Tokyo. 



A Conventions 

Since we heavily rely on the U(oo) representation theory and symmetric polynomials, let 
us explain briefly the notation following [26]. Recall that, via the Frobenius formula, the 
Schur polynomial labeled by a Young tableau R is expressed as (for an N x N unitary 
matrix U) 

/pirn t ii \^ Xn(C(k)) det U j i+N ~ l 

(R\U) = TrU = ^ T i (U) = —— N=r , (A.l) 

R t z k detu j 



9 As pointed out in [37], i and m here assign an U(N) representation via a rectangular Young tableau 
of m rows and I columns. 
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where u^s are eigenvalues of U, while ki is the i-th row length of R. In addition, xr(C(&)) 
is the character w.r.t. R evaluated at the conjugacy class C(k) of the symmetric group 
Si where 

/ = £j%, k = (ki,k 2 ,---) (A.2) 

3 

and 

oo 

?m = ikw)**, z % = n%! a (a.3) 

The character i s referred to as (R\k) if one quotes the 2d CFT techniques in 

(1231) and (l2~T0l) to define what follow: 



1$ = IK«-i)*'l°>> l^> =^)Il6- mi -jc_ n4 _i|0>, 

j=l i=l 

J2-Jk)(k\ = 1, (R\R') = 5r, ri , e(R) = ex P tn(j2^ + \r(r-l)), 

u k i=l — 



(A.4) 



where r is the number of diagonal boxes of R, while integers m; and ni (i = 1, • • • , r) are 

n^i = hi — i, rii = Vi — i. (A. 5) 

Here, /i, (uj) labels the i-th row (column) length. These tools altogether facilitate another 
expression for a Schur polynomial, i.e. 

s R (u) =TrW= (R\u), \u) = exp £^a_ n |0>, w n = TrW". (A.6) 

R n=l 71 

B Non-abelian instanton counting 

For M = 2 S77 (m) gauge theory, the non-perturbative part of Nekrasov's formula is 
Z *« = '£aM J] o,s-a r + h(k Sti -/;V. .!-')_ \R\ = Y,\R i \, (B.l) 

£ (s,i)^(r,j) a * a r + h U I) i=l 

where A is the dynamical scale and R consists of m Young tableaux located at a r 's|ll|. 
According to [38], given an arbitrary Young tableau A, we can divide it into m subparti- 
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tions C R and convert its Plancherel measure into fIB.ll) as follows: 



1 - 2 - lim ]\( ''' : h ' J ' J h 

m Ns h — h ■ A- n 7 m ^ iVr TIA _ TXr _|_ L . _ L ._)_„•_„• 

TT TT ,4 > J ■* TT TT TT ■* 

i y y y y w a - w r + 3 - % 

s=l i^tj J s=£r i=l j=l 4 i • J /g 



a s — a r + ft(j — i) 



where a r 's denote vevs of complex scalars in vector multiplets of unbroken U(l) m ~ l C 
SU(m). 
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